Ion specificity and the theory of stability of colloidal suspensions 
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A theory is presented which allow us to accurately calculate the critical coagulation concentration 
(CCC) of hydrophobic colloidal suspensions. For positively charged particles the CCC's follow the 
Hofmeister (lyotropic) series. For negatively charged particles the series is reversed. We find that 
strongly polarizable chaotropic anions are driven towards the colloidal surface by electrostatic and 
hydrophobic forces. Within approximately one ionic radius from the surface, the chaotropic anions 
loose part of their hydration sheath and become strongly adsorbed. The kosmotropic anions, on 
the other hand, are repelled from the hydrophobic surface. The theory is quantitatively accurate 
without any adjustable parameters. We speculate that the same mechanism is responsible for the 
Hofmeister series that governs stability of protein solutions. 
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All of biology is specific. The ion channels which con- 
trol the electrolyte concentration inside living cells are 
specific to the ions which they allow to pass. The tertiary 
structure of proteins is sensitive to both the pH and to 
specific ions inside the solution. It has been known for 
over a hundred years that while some ions stabilize pro- 
tein solutions, often denaturing them in the process, oth- 
ers lead to protein precipitation. In the fields as diverse as 
biophysics, biochemistry, electrochemistry, and colloidal 
science, ionic specificity has been known — and puzzled 
over — for a very long time. It has become known as 
the "Hofmeister effect" or the "lyotropic" series of elec- 
trolytes, depending on the field of science. The tradi- 
tional physical theories of electrolytes completely fail to 
account for the ion specificity. The Debye-Hiickel (DH) 
theory of electrolytes and the Onsager-Samaras theory of 
surface tensions treat ions as hard spheres with a point 
charge located at the center [l[ . The cornerstone of col- 
loidal science, the Derjaguin, Landau, Verwey and Over- 
beek (DLVO) theory of stability of lyophobic colloidal 
suspensions is based on an even simpler picture of hard 
sphere-like colloidal particles interacting with the point- 
like ions through Coulomb potential. The DLVO theory 
showed that the primary minimum of colloid-colloid in- 
teraction potential — arising from the mutual van der 
Waals (dispersion) attraction — is not accessible at low 
electrolyte concentrations because of a large energy bar- 
rier. When the electrolyte concentration is raised above 
the critical coagulation concentration (CCC), the barrier 
height drops down to zero leading to colloidal floccula- 
tion and precipitation. The DLVO theory, however, pre- 
dicts that the CCC should be the same for all monovalent 
electrolytes, which is clearly not the case @-[f|. In fact 
it has been known for a long time that the effectiveness 
of electrolyte at precipitating hydrophobic colloids fol- 
lows the lyotropic series. For positively charged particles 
the CCC concentration of sodium thiocyanide is an or- 
der of magnitude lower than the CCC of sodium fluoride. 



Even more dramatic is the variation of the CCC with 
the sign of colloidal charge. The CCC's predicted by the 
DLVO theory are completely invariant under the colloidal 
charge reversal. Therefore, changing the charge from say 
+0.04 to -0.04 C/m 2 , results in exactly the same CCC. 
This is completely contradicted by the experiments which 
find that under the same charge reversal, the CCC's can 
change by as much as two orders of magnitude for ex- 
actly the same electrolyte. Similar dramatic effects are 
observed for protein solutions for which some electrolytes 
are found to "salt-in" while other "salt-out" the same pro- 
tein. In view of the fundamental importance of ion speci- 
ficity across so many different disciplines, there has been 
a g reat effort to understand its physical mechanisms @- 
14[. A successful theory should be able to quantitative 
predict the CCC's of colloidal suspensions and shed new 
light on the specific ion effects so fundamental to modern 
biology. Construction of such a theory is the subject of 
the present Letter. 

Over the last twenty years there has been a growing 
realization that ionic polarizability — the effective rigid- 
ity of the electronic charge distribution — plays an im- 
portant role in coding the ionic specificity 15 -Tj^. The 
work on surface tension of electrolyte solutions showed 
that near an air-water interface ions can be divided into 
two classes: kosmotropes and chaotropes [lflL l20|. The 
kosmotropes remain strongly hydrated, and are repelled 
from the interface. On the other hand, chaotropes loose 
their hydration sheath and redistribute their electronic 
charge so that it remains mostly hydrated. This way 
chaotropes gain hydrophobic free energy at a small price 
in electrostatic self-energy. For hard non-polarizable ions 
of the Debye-Huckel-Onsager-Samaras theory, the hy- 
drophobic forces are just too weak to overcome the elec- 
trostatic self-energy penalty of exposing the ionic charge 
to the low dielectric air environment, forcing these ions 
to always stay in the bulk, contrary to simulations and 
experiments. 
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It is natural to ask if the same mechanism is also re- 
sponsible for the ionic specificity of the CCCs and for 
the stability of protein solutions. The similarity between 
the air-water interface and a hydrophobe- water interface, 
however, is not so straight forward [2l|. The simulations 
find that contrary to the air- water interface, the mean 
water density is actually larger near a hydrophobic sur- 
face than in the bulk [22j. On the other hand, the in- 
terfacial water is "softer" than bulk water, so that the 
hydration layer is very compressible 22]. In particular, 



the simulations find that the energy cost of placing an 
ideal cavity of radius a near a hydrophobic surface is half 
that of having it in the bulk [23]. For small cavities, 
the bulk energy cost scales with the volume 21, [23| as 
Ubuik = va 3 , where v = 0.3fcsT [24j . Moving an ion to 
the surface, therefore, gains U cav — — § fl3 of hydrophobic 
energy. For kosmotropic ions this energy is too small to 
compensate the loss of the hydration sheath and for ex- 
posing part of the ionic charge to the low dielectric envi- 
ronment of the surface. For highly polarizable chaotropic 
ions the cavitational energy gain compensates the elec- 
trostatic energy cost, since these ions can shift most of 
their electronic charge to still remain hydrated. Further- 
more, if there are no water-molecules between an ion and 
a surface, dispersion interactions — resulting from quan- 
tum electromagnetic field fluctuations — come into play. 
Dispersion forces are strongly attenuated by water. For 
example, the Hamaker constant for polystyrene in wa- 
ter is seven times lower than it is in vacuum 25]. The 



London-Lifshitz theory predicts that the dispersion inter- 
action between an ion and a surface decays with the third 
power of separation and is proportional to the ionic po- 
larizability. Although there is no exponential screening 
of dispersion force by electrolyte, in practice, it is very 
short ranged and is only relevant for chaotropic ions in 
almost a physical contact with a hydrophobic surface, see 
Fig.ffl 
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FIG. 1: Illustrative representation: while the chaotropic ions 
loose part of their hydration sheath near a hydrophobic sur- 
face and become adsorbed, the kosmotropic ions are repelled 
from the surface by the hydration layer. 

We shall work in the grand-canonical ensemble. The 
colloidal particles will be treated as planar surfaces of 
surface charge density a with the electrolyte in between. 
To account for the curvature we will use the Derjaguin 
approximation 25] . The static dielectric constant of the 



medium (water) is e w . The Bjerrum length is defined 
as Xb — (3q 2 /e w , where (3 = l/fc^T, and is 7.2 A, in 
water at room temperature. The system — two plates 
with electrolyte in between — is in contact with a salt 
reservoir at concentration ,05. The number of cations 
and anions N± per unit area in between the surfaces, 
as well as their spatial distribution will be calculated by 
minimizing the grand-potential function. Note, that we 
do not require the charge neutrality inside the system — 
the external electrolyte will screen the excess charge. Of 
course, if the system is in vacuum, the charge neutrality 
will be enforced by the formalism developed below. 
The grand-potential function is 
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where A is the surface area, x is the distance measured 
from midplane, L is the separation between the plates, 
p±(x) are the ionic concentrations, A = qN- — qN + — 2a 
is the deviation from the charge neutrality inside the 
system, and U±(x) are the interaction potentials be- 
tween the ions and the surfaces. The first two terms 
of eq. <m are the entropic free energy, the third term 
is the electrostatic energy inside the system, the sixth 
term is the electrostatic penalty for violating the charge 
neutrality, and the last term accounts for the work that 
must be done against the pressure of reservoir. The 
electrostatic penalty (sixth term) is calculated using the 
DH equation (linearized Poisson-Boltzmann equation) 
and accounts for the screening of the excess charge by 
the external reservoir. The inverse (screening) Debye 
length of the reservoir is k w — ^/8ttXbPs- The equi- 
librium values of N± and the density profiles p±(x) 
are determined by numerically minimizing the grand- 
potential function. Using the Gauss law, the electric field 

is E{x) = — - Jq {p + (x') — p_(x')] dx' . Minimizing the 

grand-potential function for a fixed N±, we find 



At 0*0 
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where the electrostatic potential is (f>[x) = — f Q dx'E(x'), 
O is the Heaviside step function, a± are the hydrated or 
unhydrated (bare) ionic radius, depending if the ion is a 
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kosmotrope or a chaotrope. Substituting the ionic den- 
sities into Gauss' law, yields an integral equation for the 
electric field which is then solved numerically. We pro- 
ceed as follows: for a trial number of cations and anions, 
N±, we calculate the electric field and the ionic density 
profile. Once these are obtained, we substitute them back 
into the free energy functional ([J). The trial values of 
N± are varied until the minimum of the grand-potential 
is found. In practice, we couple the integral equation 
solver to a minimization subroutine. 

The van der Waals interaction between two surfaces, 
is given by [Hj]: H(L) = -A™ p /l2irL 2 , where A™ p is 
the Hamaker constant for polystyrene in water j25[. To 
account for the finite radius of colloidal particles, we use 
the Derjaguin approximation [25j ]. The total interaction 
potential between the two latex spheres of radius R c is 
then U tot (L) = ttR c dl [9.(1) - fi(oo) + H(l)\. 




L(A) 



FIG. 2: Particle-particle interaction potentials in units of 
ksT/A, for suspensions with sodium-halide salts. The col- 
loidal surface charge is 0.04 C/m 2 . The electrolyte concen- 
tration is ps = 20 mM. While suspensions containing NaF, 
NaCl, and NaBr at this concentration are stable - there is 
an energy barrier for the primary minimum - suspension with 
Nal is unstable. 

We begin by studying the colloidal suspensions with 
sodium salts of kosmotropic anions, IO^ , F~, BrO^ , 
Cl~ . For these salts both cations and anions are strongly 
hydrated and are repelled from the colloidal surface. 
Since the static dielectric constant of latex is less than of 
water, in addition to the hard-core repulsion at contact, 
the ions also feel a repulsive charge-image interaction. 
To explicitly calculate this potential one needs to resum 
an infinite number of images. This can be done using 
the theory developed in Ref. [26[, yielding the ion-image 
interaction potential Wi m (x) . 

For kosmotropic ions the latex-ion interaction potential 
has a particularly simple form — hard core repulsion at 
one hydrated ionic radius from each surface plus the ion- 
image interaction. The hydrated radii were taken from 
Nightingale [13] and are the same as used in our previous 
work on surface tensions 0, [20| . The interaction poten- 
tial between the two colloidal particles, Utot, can now be 
calculated, Fig. [21 The potential has a primary minimum 



at L = 0, followed by an energy barrier. The CCC is de- 
fined as the concentration ps for which the barrier goes 
down to zero. In Table U we present the calculated val- 
ues of CCC for positively charged colloidal suspensions 
and compare the theoretical results with the experimen- 
tal measurements. A very good agreement between the 
theory and experiments is found without any adjustable 
parameters. To further test the accuracy of the theory we 
have calculated the CCC for negatively charged colloidal 
particles with a = —0.061 C/m 2 in a suspension with 
NaCl [3]. In this case the experimental value of CCC 
was found to be 140 mM, while we obtain 134.2 mM. 

For chaotropic anions the situation is significantly 
more complicated. As mentioned earlier these ions can 
shed their hydration sheath and adsorb to the colloidal 
surface, gaining hydrophobic free energy. In addition, 
diminished hydration leads to strong dispersion interac- 
tion between latex and a chaotropic anion. The London- 
Lifshitz theory predicts that the ion-particle dispersion 
potential has the form U d i S (x) ps B[ ^ L ^ +x ^ + {L /2_ X )A , 
where B is the constant related to ionic excess polar- 
izability and the ionization potential, as well as the di- 
electric properties of polystyrene and water. Unfortu- 
nately these quantities are not known. However, we 
can get a reasonable approximation for B using a sim- 
plified Hamaker theory (25j . Within this approxima- 
tion we find B — 2A c sa/%, where A e g is the effective 
Hamaker constant for a metal-polystyrene interaction 



and a is the ionic polarizability 18[. In vacuum, us- 
ing a standard relation for Hamaker constants, we obtain 
A v e g = y A pp A mm , where A 

mm are the constants 

for metal-metal and polystyrene-polystyrene interaction, 
resulting in A v cii = 17.776 x 10~ 20 J. On the other hand if 



the ion is fully hydrated A™ s 



. A w where A w 

pp mm, v ^ mm 

and A pp are the Hamaker constants for metal-metal and 
polystyrene-polystyrene interaction in water [2oT |. yield- 
20 J. A chaotropic ion near a latex 



ing A™ s = 6.245 x 10" 
surface, however, is only partially hydrated, so that the 
effective Hamaker constant should be intermediate be- 
tween Alff and A™ s , which we take to be the arithmetic 
average, A cS = [A™ s + A v cS ) /2. 

As the ion moves away from the surface, it creates 
an additional cavity from which water molecules are ex- 
cluded. This, once again carries a hydrophobic free 
energy cost. Since the volume of this new cavity is 
small we can, once again, use the volumetric scaling of 
the cavitational free energy to estimate its cost. This 
solvophobic free energy will grow linearly as /3U so i(x) = 
—va 3 _/2 + (3/4)b i a 2 _(L — x — a-) up to the maximum sep- 
aration 5a_/3 from the surface, after which distance it 
will be zero — the cavitational energy will take its bulk 
value. 

Thus for chaotropic ions, in addition to the ion-image 
interaction, the solvophobic and the dispersion energies 
must be taken into account: U-{x) = Wi m (x;a—) + 
U so i{x) + Udis(x). Finally, unlike the kosmotropes - 
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which are limited by their hydration — the chaotropes 
can come to the colloidal surface up to one bare ionic 
radius. The ionic radii and polarizabilities for I~ , NO^ , 
Br~, ClO^ and CZ07 are the same as used in our pre- 

3 



The CCC's are 



vious work on surface tensions 
once again obtained by requiring zero energy barrier for 
the primary minimum. Table U compares the theoret- 
ical predictions with the experiment. The good agree- 
ment in this case, however, might be fortuitous consid- 
ering the crudeness of our treatment of the dispersion 
interaction. In order to test that this is not the case, 
we have calculated the CCC for a sodium-nitrate salt 
in suspension of negatively charged colloidal particles of 
fj = —0.061 C/m 2 , the CCC of which is known experi- 
mentally to be 170 mM Q. The present theory predicts 
163 mM, suggesting that the good agreement between 
the theory and experiment is not a coincidence. 



TABLE I: Comparison between the calculated CCC's and the 
experimental values. The colloidal surface charge density is 
+0.04 C/m 2 from reference [5|. Judging from the spread of 
experimental data, the experimental error is approximately 
±10 mM. 



Ions 


theory (mM) 


experiments (mM) 


JO3- 


79.5 


85 [5] 


F~ 


78 


80 [5] 




71.7 


* 


cr 


70 


70 [5], 90 [3] 




31 


36 [5], 35 [3] 


ClO^ 


25.08 


* 


Br" 


24.3 


45 [5] 


cioi 


20.9 


* 


r 


14.8 


18 [5] 


SCN~ 


8.625 


12 [5], 27 [3] 



The case of thiocyanide SCN~, is somewhat more 
complicated since this ion can not be modeled as a 
sphere, but is rather a cylinder of radius 1.42 A and 
length 4.77A (2^| . The transverse polarizability of SCN~ 
is taken to be 3. OA 3 (29J. It is energetically favorable 
for this ion to be adsorbed flat onto colloidal surface. 
The cavitational energy can be calculated similarly as 
above. With these parameters we find that SCN" is very 
strongly adsorbed at the colloidal surface, see Table fl] 
Finally, we note that for negatively charged hydropho- 
bic surfaces, the CCC's follow the reversed Hofmeister 
series. For example, for a suspension of particles with 
a = -0.04 C/m 2 the CCC's range from 70 mM for JO3 
to 239 mM for SCN~. Unfortunately no experimental 
data is available in this case. 

We have presented a theory which allows us to 
quantitatively predict critical coagulation concentrations. 
For positively charged particles, the CCC's follow the 
Hofmeister series. The only significant deviation is Br~ . 



If the experimental data is correct, this would suggest 
that the values of polarizability for halogen ions quoted 
in the literature are overestimated. This seems to be 
consistent with the recent conclusions of ah initio simula- 
tions [30j. For negatively charged particles the Hofmeis- 
ter series is reversed. In view of the success of the present 
theory it is now reasonable to hope that a fully quanti- 
tative understanding of the Hofmeister effect for protein 
solutions might be in sight. 

This work was partially supported by the CNPq, 
INCT-FCx, and by the US-AFOSR under the grant 
FA9550-09-1-0283. 
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